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Abstract 

Using the Karpman-Solov'ev method we derive the equations for the 2- 
sohton adiabatic interaction for solitons of the modified nonUnear Schrodin- 
ger equation (MNSE). Then we generalize these equations to the case of 
A'^ interacting solitons with almost equal velocities and widths. On the 
basis of this result we prove that the N MNSE-soliton train interaction 
(A'^ > 2) can be modeled by the completely integrable complex Toda chain 
(CTC). This is an argument in favor of universality of the complex Toda 
chain which was previously shown to model the soliton train interaction 
for nonlinear Schrodinger solitons. The integrability of the CTC is used 
to describe all possible dynamical regimes of the A''-soliton trains which 
include asymptotically free propagation of all A*' solitons, A^'-soliton bound 
states, various mixed regimes, etc. It allows also to describe analytically 
the manifolds in the 4A'^-dimensional space of initial soliton parameters 
which are responsible for each of the regimes mentioned above. We com- 
pare the results of the CTC model with the numerical solutions of the 
MNSE for 2 and 3-sohton interactions and find a very good agreement. 

1 Introduction 

The analytical description of the dynamics of picosecond solitons in single-mode 
nonlinear fibers is based on the nonlinear Schrodinger equation (NSE) |^. The 
NSE serves as a very good integrable model admitting comprehensive investiga- 
tion in the framework of the inverse spectral transform (1ST) 1ST provides 
the complete analytical description of the soliton interaction in a generic case of 
asymptotically free N solitons moving with pair- wise different velocities ^ . 
On the other hand, the practically important case, especially in a soliton-based 
fiber transmission, deals with the so-called A^-soliton trains, i.e., with an or- 
dered sequence of N {N > 2) solitons which are spaced apart almost equally 
and have almost (or exactly) equal amplitudes and velocities. In a number of 
recent papers ^, |^, ^ , an effective formalism was developed for studying the 
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dynamics of well-separated NSE solitons within the A'^-soliton train. This ap- 
proach is based on a generalization of the two-soliton quasiparticle method by 
Karpman and Solov'ev |^ to the case of N solitons. In the framework of this 
approach, the soliton interaction is governed by a dynamical system for AN soli- 
ton parameters. Such an approximation is called adiabatic because interaction 
between the solitons is displayed as a slow deformation of their parameters, a 
possible presence of radiation being ignored. It is important to realize that the 
above generalization from two to iV solitons is nontrivial because of lack of the 
superposition principle for the nonlinear dynamical system. 

Under some additional restrictions imposed on the soliton parameters, which 
ensure the validity of the adiabatic approximation, the above dynamical system 
is reduced to the complex Toda chain (CTC) equations with N nodes ||ic|] . 
Extensive use of the fact that the CTC is a completely integrable model permits 
to classify soliton parameter regions with different asymptotic regimes of the N- 
soliton train H, ^, 0, || . It was also shown in j|] that the CTC can be associated 
with any equation from the NSE hierarchy. 

One of the purposes in the optical fiber soliton communication is to achieve 
as high of a bit rate as possible. A natural way in this direction is the use of 
shorter optical pulses. It should be noted, however, that when dealing with 
ultrashort optical pulses with duration < 100 fs, the NSE should be modified to 
take into account some additional effects, such as the nonlinearity dispersion, 
the intrapulse Raman scattering and the higher-order dispersion [0. As a rule, 
the extra terms added to the NSE violate its integrability. On the other hand, 
if these additional terms are small, the IST-based soliton perturbation theory 
is usualy treated as the relevant method to account for their influence on the 
soliton behavior 0, |l2|, |l3|. 

It is remarkable that adding a term accounting for the nonlinearity dispersion 
to the NSE preserves the integrability of the equation. In other words, the 
modified nonlinear Schrodinger equation (MNSE) 



is still integrable by means of 1ST, though the associated spectral problem (the 
so called Wadati-Konno-Ichikawa spectral problem ||lj], or quadratic bundle) 
does not belong to the familiar Zakharov-Shabat class. The parameter a in 
(0) governs the strength of the nonlinearity dispersion. The case a = corre- 
sponds to NSE. Thereby, the effect of the nonlinearity dispersion is considered 
non-perturbatively in (|l|). Moreover, we have to stress that it is the completely 
integrable model (|l|) that should be considered as a true starting point for an- 
alytical investigation of subpicosecond soliton dynamics. Indeed, it was shown 
in (Tsl that numerical simulation of the soliton propagation according to the 
MNSE (|l|) revealed various kinds of dynamical behavior which cannot be ac- 
counted for by treating the nonlinearity dispersion term of the MNSE (]l|) as a 
perturbation term in the NSE. Analogous idea in treating the perturbed NSE 
was developed by Kodama and Hasegawa in ||l^ . There the NSE with pertur- 
bations like the third order dispersion, nonlinear gain and nonlinear dispersion 
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was treated as a perturbed higher order NSE. 

The relevance of Eq. M)to the problem of ultrashort pulse propagation in 
fibers was demonstrated in . MNSE is also used in plasma physics 

and is relevant for description of a deformed continuous Heisenberg ferromag- 
net . It is the Alfven waves in magnetized plasma where the first successful 
application of 1ST to the quadratic bundle was achieved on an example of the 
derivative NSE which is Eq. (|l|) without the last term. Both equations 
are interrelated by a gauge-like transformation, see, for example, |^ , |2^ . The 
soliton solutions and the Hamiltonian structures of the MNSE were obtained 
for the first time in A^-soliton solutions were further rederived by differ- 

ent methods: by 1ST using the above relation with the derivative NSE ||2^ , 
by Backlund and Darboux transformations |2^, by technique of determinant 
calculations [Q, by the Hirota method |27 , by the 9- method |2^. It should 
be noted that the solutions obtained in these papers refer to the general case of 
asymptotically free solitons and being exact were too complicated for practical 
use. 

Recently, a novel parametrization for the MNSE solitons was proposed within 
the framework of the Riemann-Hilbert formulation of 1ST . The convenient 
parametrization of the MNSE soliton facilitated the development of an effective 
adiabatic soliton perturbation theory for the MNSE which is able to take into 
account non-zero terms in the right-hand side of (|l|), see p9[ . 

The next natural step is to derive dynamical equations of the Karpman- 
Solov'ev type for the adiabatic evolution of the soliton parameters for the MNSE 
iV-soliton train. Several questions arise in the process of solving this problem. Is 
it possible to associate an Ai"-node chain model, like the CTC, with this dynam- 
ical system? Will this chain model be different from the CTC and, therefore, is 
the CTC valid only for the NSE hierarchy? How well do the numerical simula- 
tions of the MNSE with adiabatic A^-soliton train initial conditions agree with 
the chain-like model predictions? All these questions will be answered below. 

The purpose of this paper is to derive a dynamical system for the AN soliton 
parameters for the MNSE 7V-soliton train. To this end we will generalize to the 
quadratic bundle the similar investigations performed for the NSE. In the next 
section we apply the Karpman-Solov'ev approach to the MNSE (|l|) and derive 
the dynamical system for the 2-soliton train. In Sec. 3 we show how this result 
can be generalized to the MNSE A^-soliton train for N > 2. We show in Sec. 
4 that after some additional assumptions the corresponding dynamical system 
for the soliton parameters acquires the form of the CTC. Thus we find that the 
CTC is characteristic not only for the NSE hierarchy ||], but has wider field of 
applications. This is an argument in favor of its universality. 

In Section 5 we show how the integrability of the CTC can be used to 
determine the dynamical regimes of the A^-soliton trains. We demonstrate on 
the examples of = 2 and A^ = 3 how one can describe analytically the 
manifolds in the 4Af-dimensional space of initial soliton parameters which are 
responsible for the: i) A^-soliton bound state regime; ii) asymptotically free 
regimes; iii) various mixed regimes, etc. Although the analysis follows closely 
the ideas developed in ^, Q the description of the corresponding manifolds 
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differs from the ones for the NSE soliton trains. The reason for this Hes in the 
fact that the CTC fields Qj{t) are parametrized in a different way; in particular, 
Im Qj depend not only on the soliton phases 5j (as is the case for the NSE) but 
also on the soliton amplitudes. 

In Secton 6 predictions of the CTC model are compared with the numerical 
results from the MNSE and find an excellent match for most regimes with N ~ 2 
and = 3. We found some disagreement between the CTC and numerical 
MNSE solution in the regimes when CTC predicts a very slow soliton separation. 



2 2-soliton interactions for the MNSE 

First of all we summarize the basic results concerning the soliton solution of the 
MNSE (|l|) 1^ . This equation admits the Lax representation 

= ~-(k^-]][a3,^] + 2ikQ<i>, (2) 



a \ 4^ 

Ai i \ 

— k^Q + 2ik'^Q'^(j3 kQ + kQ^a^i - 2iakQ^ 

a a J 

Here the Hermitian matrix Q — ^ ^ stands for the potential of the spectral 

problem A; is a spectral parameter. There exist various paramctrizations of 
the soliton solution of the MNSE, the first one having been proposed in | |2^ . 
We follow here the parametrization given in ||2^ which was proven to be useful 
for practical calculations and admits a simple (though nontrivial) reduction to 
the NSE for a ^ 0. Namely, the MNSE soliton solution is written as follows: 

"^^^'*^"*a(fce^ + fce-^)2' " ^"^^ 
Here k ^ ku — ikj, fc/ > 0, A = 4fc^ — 1 = /i — 

a V 
m^~^t + ^^o): S{t)^^{^,' + iy^)t + d^o)- (4) 

It should be stressed that the soliton (^) is not of the hyperbolic-secant type 
with a real argument, characteristic for the NSE. It is specified by four real 
parameters /i, v, ^(g): ^ind (5(o) with {—fi/a) being the soliton velocity, {a/u) is 
its width, (5(0) and ^(q) are initial phase and position of the soliton. To carry out 
the limit reduction to the NSE, one should decompose the spectral parameter 
in the following manner: 

fc=^-f + «^0, (5) 
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which gives {-li/a) 2^i^^^ and (a/i^) {2p^^^)-'^, as should be. 

If there is a smaU perturbation in a system described by the MNSE, we wiU 
deal with a perturbed MNSE: 

iut + ^Uxx + ia + \u\^u = r{x, t), (6) 

where r(x, t) describes a functional form of the perturbation. In what follows we 
will restrict ourselves to the adiabatic approximation of the soliton perturbation 
theory. In other words, we suppose that a perturbation causes a slow variation 
of the soliton parameters only. The evolution equations for the perturbation- 
induced soliton parameters are given in [p9[ . Here we write them in terms of 
the parameters (^. The key equation has a very simple form: 

dt 2 {ke-' + ke^Y ' ^ ' 

where R± = exjp[~i<j){z,t)\r{z,t) ± exp[i(/)(— z, t)]f(— z, t). Taking into account 
we obtain 

(fee- + ke^-Y + ^ ' 

Evolution of ^ and 5 is given by the following formulae: 

t 



at a av \jn J at 



r°° dzi?_ 



z (fc'^e^ + t^e-"") + — ike'' + ke" 



,(10) 



' / d<'/i(i') /i— - 



dt 2oP- a^v \Jq ^ V V di dt 

\k^\ (fee ^ + fce^ 



1 /T , , . T ,s jz 



(Afce^ - A/ce"^) + — Afc-^e^ + Afc-^e" 



(11) 



It should be noted that for the symmetric perturbations obeying the condi- 
tion exp[— i)]r(z, i) = exp[z0(— z, t)]f(— z, t\ i.e., i?_ = 0, the complicated 
integrals in the right hand sides of (|l^) and (^ij) disappear. 

Now we have all the necessary information to derive the Karpman-Solov'ev- 
like dynamical system of equations for the adiabatic interaction of two well- 
separated MNSE solitons. Below we will formulate more precisely the condition 
of sufficient separability of solitons. We suppose that a two-soliton solution to 
the MNSE (§) is well approximated by the sum of two MNSE solitons: 

u{x,t) = Ui(zi,t) + U2{z2,t), (12) 
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where Uj{zj,t), j = 1,2, is given by (1) with 



where we took into account for the possible evolution of and Vj . Substituting 
now (|l2|) into the MNSE (|l|), it is easy to see that, because of the nonlinearity, 
each soliton feels the presence of the other one and the interaction is described 
by the perturbed MNSE 

iujt + ]^uj^^ + ia {\uj\^Uj)^ + |ujf uj = rj, (13) 

where 

Tj = -ia (2|ujpM3_j + WjW3-j)^ - {2\uj\^U3-j + v^jU^^j) ■ (14) 

It should be stressed that the perturbation (^^ arises effectively as a result of 
treating two-soliton solution as the sum ( |l2[ ) of the one-soliton ones. 

Now we formulate the conditions which provide the representation ( p^ ) as a 
two-soliton solution of the MNSE (|^) . At first we express Z2 in terms of zi , 

^2 = IH zi + — (^2 - Ci)- 

V J a 

We suppose that sohtons have almost equal widths, i.e., 

« 1, (15) 

where vq — (l/2)(i^i + ^2)- Hence, we have 

Z2- 21^^(6 -6)- (16) 
a 

Calculation of the overlap integral | 1/1(21, t)M2(z2, i)dx| (or, equivalently, 

da;|Mi(zi, t)u2{z2, t)\) gives an expression containing the factor exp[— (^'o/Q!)(^2- 
^l )] = e for ^2 > Ci ■ Just this exponential factor determines a measure of over- 
lapping neighboring solitons. We take in the following 

-I6-Cil»i (17) 

a 

(or e <C 1) which means weak overlapping between the solitons. 

Let us consider now the phase difference <j)2 — <t>i — {1^2/ 1^2) Z2 — {f^i/vi)zi -\- 
82 — 5i. Accounting for ([l5| ) and (|l^) we may write 



V2 



At2 - ( IH — ) Ml 



zi + — — (6 - 6) + 52- 5i- 
a V2 
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Since we consider solitons moving with small relative velocities we assume: 

« 1. (18) 

I'D 

Then the phase difference will not contain the z-dependence. Furthermore, 

a V2 a \ V2 ) 

As the last condition we suppose 

(19) 

hence, the phase difference takes the form 

J, A. ^0 , 



cs. 



-(6-a) + '52-<5i. 



Therefore, the conditions (^5|), ([T^), (|T^ ), ( [T9| ) provide a possibility to consider 
a two-soliton solution of the MNSE (pj) as the sum of the form (|l|). 

To derive the Karpman-Solov'ev-like equations for the soliton parameters, 
we use Eqs. ([7|)-(pl|) with the perturbation (p^). Accounting for the above 
conditions, we obtain after simple but tedious calculations (j — 1,2): 



^ = (_l).4fc,. !^^3-^^_A-.^ (20) 

dt ^ ' ^ \ a I ks-j ^ ' 



-Mj + (-1)^— i'j^^a-j / dtVj(0 



dt a 



%3-, [(1 + -f) (1 - 2w^) + 4^s,] ^«>3-, ie-^ 

[(1 + w^) (1 - 2w^) - 4is,] wp,^,-^eA e-^, (23) 

K3-J / 



dJ, 

dt 2a 
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^j'^s-j [(1 + ^5 ■) (1 - 2i5,') + 4zs,] u;>3-, 



[(1 + (1 - 2w^) - 4z.,] w^^,_^-^eAe-^. (24) 



Above we used the notations 



kj a 



V' = — (6 - Ci) + (52 - <5i, Sj = - arctan 



a 2 1 + /ij 

where the last relation follows from \j = 4fc| — 1 — fij — ivj. 

Equations (^-(p4[) are the analog of the Karpman-Solov'ev equations in 
the case of the adiabatic interaction of two well-separated MNSE solitons and 
reduce to the NSE dynamical system in the limit (||). 

The dynamical system (pl|)-(p^) is rather complicated and needs further 
simplification to perform its analytical investigation. Integrable approximation 
is of special importance, and finding such an approximation is one of our pur- 
poses. But first of all we will generalize these equations to the case of iV MNSE 
solitons. 



3 A'^-soliton train interactions for the MNSE 

Since the Karpman-Solov'ev- like dynamical system is nonlinear, it does not 
allow the superposition principle. It is physically clear because in the case 
of the A^-soliton train with > 3 a middle soliton will be influenced by its 
neighbors from both sides. Hence, it is not possible to describe the interaction 



of A > 3 solitons within the framework of two-soliton interaction like (21 )-(|24|). 

The first remark we should keep in mind is that the interaction force between 
the solitons is of the order of their overlap. Therefore, we can take into account 
only the nearest neighbor interaction. Indeed, for the A-soliton train we assume 
that 

N 

i=i 

where Uj is the MNSE soliton (||) whose center of mass is located at . Assume 
that £,1 < £,2 < ■ ■ ■ < £,N ■ Inserting this ansatz into the cubic term of the MNSE 
(1) gives 



N 

UjUlUn- (26) 

j = l j^l i=^l=^n 
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Straightforward analysis shows that the integrals in ([7|)-(pl]) corresponding to 
each type of terms in (26) are of the following order of magnitude: 



UjUiUn, j<l<n ^ C'(e"l-'^'l"l'""l). 

Here e is of the order of exp[— (1/0/0)1^^— for — — 1. Because we keep only 
terms of the order of e, we see that it is enough to take into account only terms 
with \j — /| = 1. In other words, the 'triple' terms like UjUiUn can be neglected. 
Quite analogous is the situation with the cubic terms containing a;-derivative. 
Secondly, as in Sec. 2, we pose the conditions 

where vq = N^^ SjLi Mo = X^jLi Mi- They mean that we consider the 
chain-like configuration of N solitons with equal or nearly equal velocities and 
widths. Substituting the soliton solutions (H) into the perturbation 



E 

l=3±l 



\ia {2\uj\^ui + u]ui)^ + {2\uj\^ui + u]ui) (27) 
and calculating the integrals in (^, we obtain 



where 



, 1 for Z = j + A _ f^o,^ 



^i, = ^ (e, -i,) + bi~ 5, . (29) 



The corresponding formulae for /ij and i^j follow from (|2^) as real and imaginary 
parts. It is not difficult to derive also the equations for the rest two parame- 
ters and Sj generalizing those (|2^) and (^^ for the two-soliton interaction. 
Keeping in mind, however, our aim to formulate the equations for the adia- 
batic interaction of the MNSE solitons in the form tractable analytically, it is 
sufficient to represent the equations for and dj in the following form: 

^ = -^+0{eh (30) 
dr a 

Indeed, let us impose the conditions on the scattering data of the spectral prob- 
lem (U) which correspond to the adiabatic approximation. Just as for the NSE, 
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we require that the eigenvalues of the Lax operator are clustered around their 
mean value: 

1 ^ 

|A,-Aop^O(e), Ao = -^A,. 

Thus we obtain the estimates in ( ^ ) and (|3l|) which mean that we can neglect 
the perturbation-induced evolution of the parameters and dj as compared to 
their main (unperturbed) evolution. At the same time Sj and Wj characterize 
the initial conditions and it is important to take them into account in the right 
hand side of the equation (^). 



4 Derivation of the complex Toda chain model 

The next important step towards deriving a model of N MNSE-soliton inter- 
actions tractable analytically consists in a careful account for the terms of the 
order of e. First note that because the right-hand side of (|2|) is of the order of 
e, we may approximate kj by |fco|e~**^ where fco is the mean value 

N 



Thereby we neglect terms like Ivq — Vj\e and Juo — which due to (^) and 
( |l8| ) are of the higher order in e. Hence, eq. ( psl) is written as follows: 



dAo 4t'n / n 
dt a2 K'^ 



'93) 



(32) 



where 



Qj+i - Qj 



(e 



Mo , 



TT -I- — (Cj+1 - Cj) + - 6 J + 4sj+i -I- 4sj 



(33) 
(34) 



The recurrent relation (B3) can be solved for Qj with the result 



jn + + 6j + 5q + ^ 8sk + isj 



k=l 



The next step is to derive the evolution equation for Qj 
noted first of all that up to terms of the order of e 



1 ^ 

(35) 

It should be 



dJo 1 2 , 2\ 
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1 ^ 

i=i 

Then, in view of and (|3ll), we get 

^ = ^(^0 - - + + Mo + ^^o) (36) 

= ^ [-(Mj - Mo)^ - (i^j - I'o)^ - 2ii/o(Mj ~ ii^j)] 

Finally, keeping only the leading-order terms and replacing fj ~ 1 and gj ~ 1 
we find from ( |33| ) and (36) 

^ - 4 (eQ.+i-Q. - eO^-Q-i) , (37) 

dt^ V a / 

i. e., the CTC model. Hence we see that the CTC model arises naturally as the 
integrable limit of the Karpman-Solov'ev-like equations describing the adiabatic 
interaction of N MNSE solitons within the train of solitons with near velocities 
and widths. 

As we will see in the next section the interactions of the MNSE solitons are 
substantially different from the ones of the NSE. As it can be seen from ( |33|) 
the dependence of Qj on the soliton parameters is different from that for the 
NSE case. An important point here is that ImQj depends explicitly also on the 
amplitudes of the solitons through Sj — arctan (z/j/(a(l + ^j)). 



5 Dynamical regimes of the A^-soliton trains 

It is well known that the CTC is a completely integrable dynamical system. 
Most of the results concerning the CTC such as the Lax representation, the 
integrals of motion, explicit solutions etc. are direct consequences of the classical 
results by Toda and Moser ^ on the real Toda chain (RTC). However 
there is a qualitative difference between the RTC and the CTC when one tries 
to analyze the dynamical regimes of the two systems, see 0, |30|] . 
Indeed, the Lax representation of the CTC ( ^7| ) is of the form: 

^ = [L.M], (38) 

AT 

L{t) = (^J-^J.^- + «j(^J.J+i + , (39) 

N 

M{t) = ^AE,+i,, - Ej,,+i), (40) 
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where ^ 

T = Cot, ca = 2vlla^, Qj = -eyip{Qj+i ~ Qj)/2, 

bj = -idQj/dT = -Aj/4i^o, {Eij)in = SuSj,,, 

see (p6|). In fact, without loss of generaUty we can assume that tr L = 0. 
This can be achieved by subtracting CqI from L where Co = X^jLi (j/^ = 
SjLi bj/N. Note that Co is obviously an integral of motion for the CTC, i. e., 
dCo/dr = 0. 

The explicit solution to the CTC is given by 

q,{r)=q^{0) + \n^^, (41) 

where ^0 = 1, 

N 

^i(r) = 5]r^e-2a-, (42) 
fe=i 

AkiT)= J2 (ni...r,j2M/2(Zfc,;fc_i,...,;i)e-2(0i+-+«'J- (43) 
l<ii<...<ifc<Ar 

and 

N 

An = W^{N,N -1,...,2, i)e-2(Ci+...+C«)r "Q (44) 

k=l 

Here (j are the eigenvalues of the Lax matrix L, W{lk,---,li) denotes the 
Vandermonde determinant: 

Wilk,...Ji)= n (2Cs-2Cp), (45) 

s > p 

s,p e {h, . . . ,ik} 

and Tj are the first components rj = Vj^i of the eigenvectors: 

Lvj = CjVj, (46) 

normalized by 

AT 

(«„u,)^^(z/,-fe)' = l. (47) 

k=l 

Due to the fact that L is a symmetric matrix we find also 

N 

T.^] = l. (48) 

Using the explicit solution for Qj (t) we can estimate the asymptotic behavior 
of Qjir) for T —^ 00. 
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Such an analysis for the RTC, i. e., when Qj, aj, and bj are real, shows 
that: i) rj and are real-valued, ii) (j ^ C,k for j ^ k. Therefore one finds that 
for r — > cx) each 'particle' Qj moves uniformly with a velocity 2(j . Since Cj 
are pair-wise different we conclude that the only possible dynamical regime is 
the asymptotically free (AFR) one. 

The same considerations applied to the CTC lead however to a qualitatively 
different results. Indeed, now rj and Q — Kj + irjj become complex- valued and 
there are no restrictions on the eigenvalues Q. Then evaluating the limits of 
Qjir) for T — > oo we find that the asymptotic velocity of Qj is determined by 
2k j = 2Re<^j. As a result we have a much wider spectra of dynamical regimes. 
The reason for that is also in the fact that CTC can be viewed as a dynami- 
cal system of N 'complex' particles which are characterized not only by their 
positions Re Qj and velocities Vj = Re bj , but also by their phases and phase 
velocities; the latter are related to ImQj and Imbj. Physically speaking these 
'complex' particles have, just like the bright NLS solitons, an internal degree 
of freedom. This makes the interaction between the particles more complicated 
and as a result the number of the possible dynamical regimes increases substan- 
tially. 

The AFR which takes place if Kj ^ Kk for j ^ fc is just one of the options. 
Another option is ki = K2 = . • . = k^ = which corresponds to a bound state 
regime (BSR) of all N 'complex' particles (solitons) in the train. There is also 
a large class of intermediate or mixed regimes (MR) for which only several of 
the parameters equal. For example, if ki = K2 = K3 > K4 . . . > kat then 

the first three particles (solitons) will form a bound state while the rest iV — 3 
particles will be asymptotically free. 

Note that this variety of regimes exist in the generic case when the eigen- 
values Q of L are pair- wise different; so in the previous case we assume that 
Vi 7^ ^2 Vs- One may consider also degenerate regimes (when two or more 
of the eigenvalues Q become equal) and singular regimes (when one or more of 
the functions Qj{T) develop singularities for finite t). 

There is also another important consequence from the integrability of CTC. 
From the Lax representation one easily finds that the eigenvalues Cj the 
integrals of motion for the CTC, i.e., Cj ^re time independent. Therefore we 
can evaluate them, for example, at the initial moment t = using for this the 
initial values of the soliton parameters. Then, knowing Q and, more specifically, 
Kj we can predict the asymptotic regime of the corresponding A^-soliton train. 

We can also answer another question: describe the set of initial soliton pa- 
rameters for which the corresponding A^-soliton train will develop a specific 
dynamic regime. In other words, we can describe the set of initial soliton pa- 
rameters for which we will have, say, an A'^-soliton bound state regime. To 
describe the BSR all we need to do is to solve the corresponding characteristic 
equation 

det(L - CI) = 0, (49) 

and impose the condition ki — ki — . . . — kj\j — 0. Since the coefficients of ( ^9| ) 
and consequently Kj will be expressed in terms of the initial soliton parameters 
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we will have a set of nonlinear equations describing the BSR. Analogously, if we 
need to describe the AFR we must solve for Kj 7^ Kk for k ^ j. 

We will show how this can be done analytically for the simplest non-trivial 
cases with N = 2 and TV = 3. For generic values of N this can always be done 
by numeric means; one needs only to solve algebraic equation (^9|) of order N . 

Let us briefly describe the manifolds of soliton parameters responsible for 
each of the dynamical regimes for = 2 and A^ = 3. As it is clear from the 



above considerations, we have to solve the characteristic equation (49) and to 
express the eigenvalues Cj of L in terms of the soliton parameters. 

5.1 N = 2 case. 

For simplicity from now on we shall consider trains with zero initial velocities, 
fj.j{0) — 0, i. e., in the relevant moving coordinate system. The matrix Lq = 



b 



L(t = 0) = ^ " j with trL = is built from the initial soliton parameters: 
where 

''o = 6(0) -Ci(o), r = (52(o)-<5i(o)+4si+4s2, = (1/1(0) -t'o)/t^o- (50) 

Then 

Ci,2 = ±V^+^ = ±^\/yo + e-^^- (51) 

with 

A„.2 = 2e-'^«'-«/(2«), yo = (52) 

Acr,2 

Obviously if F 7^ 0, tt then ReCi,2 7^ and we will have an asymptotically free 
regime (AFR). 

If F = 0, then RcCi 2 = and wc have a bound state regime (BSR). 
If F = TT, then RcCi 2 = 0, i. e., we will have a BSR only provided 

Ml > Acr,2. (53) 

If |d| < Acr,2, both roots Ci,2 become real and we go into the AFR. 

It was already noted that the conditions F = 0, tt involve, besides the phases 
Sj, also the amplitudes of the solitons through sj. In particular for i^q — a — 
1 and Hj = we have Si + S2 — tt/A. Therefore two such MNSE solitons 
attract each other and form a bound state provided S2 — Si = tt and repulse 
each other (which leads to AFR) for S2 — Si = 0. Such behavior of the two- 
soliton interaction is quite to the contrary to that known for the NSE two-soliton 
interaction. 

The explicit solution to the CTC with A = 2 is of the form 

^ / N ^ / N , cosh(2CiCot — 71) , ri 

Qi{t) ^ -Q2(t)=ln ^, 71= In-, (54) 

2Ci r2 
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where Ci is expressed in terms of the sohton parameters (pW) and 



7i-^ln^^S±^. (55) 

Obviously for F = the solution Qi{t) 

n m 1 2cos(yoCot/2 + z7io 

Qi(.t) = In , 56 



Fo = Acr,2V2/o + 1' 710 = o / 9 . • 

2 V yo + 1-2/0 

becomes a periodic function of i = t/cq with period depending on y^. 

4n 

T,,a = — ^— . (57) 

Analogously for T — tt from (|55| ) we have 



„ ... „ , 2cosh ^Acr,2^/^/g - lcot/2 - 711) 
Qi{t) = -Q2(i) = In " , 

lAcr,2\/y0 - 1 

711 = T^ln ^-^ • (58) 

VVo - 1 - 2/0 

The solution is periodic only if > 1 smd the period is 

r2.;2 = . (59) 
CoAcr,2V2/o ~ 1 

As a conclusion, the BSR for = 2 provides periodic solutions. For F = tt, 
yo < 1 we have AFR and the solution is not periodic. 

The final remark in this subsection is that for yo ^ the solution ( p^ ) 
becomes singular and blows up periodically with period 47r/(coAcr.2)- In this 
limit we have two 'equal' solitons with amplitudes Vj = vq = I with phase 
difference tt. 

5.2 N = 3 case. 

For the case of the three-soliton train with zero initial velocities the matrix Lq 
has the form 

bi ai 
= I ai h 02 I , trLo = 

02 63 



with 
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where 

dj = , ro = ^2(0) - Ci(o) = 6(0) - 6(0): 

= <5j+i(o) - (5^(0) + 4sj+i + Asj. (60) 
Then the characteristic equation takes the form: 

C'+K + 9 = 0, (61) 

where 

16 4 

It is natural to make use of the weU known Cardano formulae for solving 
cubic equations. We first consider the cases when p and q are real. The roots 
of (|6^) are given by 

Ci =A + B, C2 = A + uj^B, C,:i^oj^A + ujB, (63) 

where 



^ - ^-f + VQ, B= - VQ, (64) 

p^ (2TTi 

If both p and q are real, then so is Q. Here we have four subcases corre- 
sponding to qualitatively different sets of roots for real p and q. 

i) Q < 0. This is possible only if p < pcr, Per = -iiq^/^Y'^. Then A = B* 
and all three roots C,j become real Qj = kj and pair-wise different: 

Ki = 2|A| cosfio, K2,3 = 2|74| cos ^fJo ± , (65) 

with flo 7^ 0, TT. Obviously this leads to AFR. If JIq = or tt then K2 — K3 and 
a MR follows. 

ii) Q > and q 7^ 0. Here both A and B are real and formula ( |6^ ) shows 
that one root Ci is real, while the other two are complex conjugate: 

ReCi = -2ReC2 = -2ReC3, or ki = -2k2 = -2k3, (66) 

which corresponds to a MR. 

iii) Q > and q = 0. Now p > 0, the cubic equation (^ij) simplifies and is 
trivially solved by 

Ci = 0, C2,3 = ±V^. (67) 
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All the roots have zero real parts which obviously corresponds to BSR. 

iv) Q = 0. If p and q are nonzero, all the rots are real and pair-wise different: 

Ci = 3g/p, C2 = Cs = -3q/2p, 

we have AFR. If p and q are zero, we get a degenerate case with all three zero 
roots. 

The symmetry in the eigenvalues leads also to a symmetry in the solutions 
of the CTC. Therefore the configuration ( |67| ) corresponds to a particular type 
of BSR's. This is due to the fact that we restricted so far both q and p to be 
real. Of course this is not necessary; moreover, from (^2|) we see that generically 
both q and p are complex. If we want to specify the soliton parameters that 
are responsible for the BSR we may also use Viette formulae which show that 
the characteristic equation ( |6l| ) will have purely imaginary roots if p is real 
and negative and q is purely imaginary. That is why we will consider also the 
configuration v) below. 

v) p = p, q = — q. In this case we have two qualitatively different possibil- 
ities depending on whether Q is positive or negative. 

Note that since q — — g we should modify our reasoning as compared to the 
above analysis. Indeed, with q ~ iq' , q' real and Q > we find that A ~ —B. 
Therefore from (|63| ) and ( |64| ) we have that all the roots Cfc satisfy Cfe = — Cfei i- 
e., are purely imaginary and BSR takes place. 

Analogously, if Q < then the roots C^k satisfy Ci = — Ci a-nd C3 = —(2 which 
leads to AFR. 

Hence, we revealed two possibilities to realize bound state regime: subcase 
iii) and subcase v) with Q > 0. 

Let us now briefly describe the sets of soliton parameters relevant to each of 
the regimes mentioned above. For definiteness we will use two configurations of 
soliton widths: 

Wl: di = -da, d2 = 0, (68) 
W2: di=d^, d2 = -2di. (69) 

The condition that p is real immediately means that 

Fi = -T2 = (70) 



Then 



p^ _l(did2-f did3 + d2d3) + ^e""°''°/"cos$, (71) 
lb z 

q = ^d,d2d, - ^e-'-^-/'^ (die'* + dae-'*) , (72) 
64 16 

<P ^ S2- (5i +4si +4s2. (73) 
Choosing the sets of widths to be Wl and W2 we get respectively: 

pW = | + |cos$, 
16 Z 

qW ^ ^^sin$, (74) 
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where eo = exp{—voro / {2a)) , and 

(2) ^ £^ + !o ^ 

^ 16 2 ' 

= _-J-__^cos$. (75) 

Case I. q = 0. The characteristic equation ( |6l| ) has the roots 

Ci = 0, C2,3 = ±V^. (76) 

From ( [T^ we get that for the Wl configuration the condition q^^^ = holds 
provided 

$ = fc7r, A: = 0,1, (77) 

which means that 

As a consequence we find that p^^-' > for fc = 0; for fc = 1 we get that p'^^ > 
only provided \di \ is greater than the critical value: 

|rfi|>A„,3, Ae,.3 = 2\/2eo. (79) 

In all these cases C2.3 are purely imaginary, i.e. these sets of parameters lead to 
BSR. 

Note that (^ means 

62 ^ Si + kn - 4si - 4,32, A: = 0,1. (80) 

If instead of we have \di\ < Acr,3 then p'^' < and the roots ^2,3 
become real. That means that taking di below the critical value we will see a 
transition from BSR to AFR. 

The same considerations applied to the W2 configuration lead to different 
results. From ( [75| ) we see that ^ holds if 

cos$ = -^, (81) 



which implies that 



and 



Mil < 260 = ^ (82) 

p(2) = ^>0. (83) 
16 

Such configurations obviously lead to BSR. If \di\ is chosen to be greater than 
the critical value in the right hand side of (|8^) we find that then q^'^'^ becomes 
purely imaginary; such situation is considered below. 
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Let us briefly treat also the case of 'equal' solitons, i.e., dj = 0. Then 



obviously q 
that if 



0, Si = S2 
TT 



S3 : 
TT 



7r/8 and p = (eo/2) cos$. As a result we find 



„ < ^' < „ 
2 2 



<S2-Si < 



then p > and we have BSR; if 



— < $ < — 
2 2 ' 



I.e. 



TT TT 

-2 <52-^i< 2, 



then p <0 and AFR follows. 

Case II. p = 0. In this case the characteristic equation (|6 



k = 0,1,2. 



(84) 

(85) 

has as roots 
(86) 



If in addition q is real then (pfl) leads to a MR; otherwise we get AFR. 
For the Wl configuration p*^^^ = means 



cos $ = — 



4 



(Acr.3)2^ 



(87) 



this is possible only if |c?i| < Acr,3. From (74) we get that q'^^^ is real and such 
configuration leads to MR. 

For the W2 configuration p^^'' = holds if 



cos ^ — 



3dl 



(Acr,3)2 ' 



which is possible only if < Ad-.s/VS. From ( [75| ) we find that q^^^ is purely 
imaginary, i.e., AFR follows. 

Case III: p = p and q = q ^ 0. This is possible only for the W2 config- 
uration, so p and q are given by ( |75|) . The resolvent of the cubic equation ( |6l| ) 
in this case takes the form: 



(p( 



2)^3 



(q(2))2 



27 



4c 

if tt: 1 r + TT 



8 12 



(89) 



where y ~ (ii/Acr,3 and c — cos$. 

It is easy to check that Q{y, c) is nonnegative for all c > —9y^/A and vanishes 
for c = and c = -~9y^/A. We have to keep in mind also that |c| < 1. Therefore 
if 9y^/4 > 1 then Q > in the whole interval —1 < c < 1. Following the 
arguments in v) above we conclude that this configurations leads to BSR. 

If we choose 

Mi|<-A„,3 (90) 



19 



then there wiU be an interval for $ (fz^), 



'/'cr < ^ < 27r — 



cr; 



(pcr = arccos — 




for which Q < 0; i.e., if holds we have AFR. 
If $ belongs to the complementary interval: 



</?cr < $ < </3, 



cr 5 



(92) 



then (5(2/, c) > and we have BSR. 

The interested reader can easily extend this studies to other relevant config- 
urations of soliton parameters. 



6 The CTC versus the numerical solutions of 



It is our aim here to compare the predictions of the CTC model with the numer- 
ical solutions of the MNSE. Since the full numerical investigation of the problem 
is a voluminous and ambitious task we limit ourselves with N = 2 and = 3 
soliton trains and fix up a = 1 and the average width vq ^ 1. 

With this choice of a = 1 the derivative term in the MNSE can not be 
treated as a perturbation to the NSE. With this choice we are able to exhibit 
the differences between the MNSE and NSE A^-soliton train interactions. As we 
mentioned above the dependence of the soliton interaction of the MNSE solitons 
on the soliton phase difference is qualitatively different from the one of the NSE 
solitons. 

Indeed, let us start with N — 2 soliton trains. The formulas from Section 5.1 
with a — 1 and 1^0 = 1 show that 'equal' solitons (i.e., solitons with equal widths) 
with phase difference 62 — Si = tt (or F = 0) attract each other. In fact this 
choice of the soliton parameters corresponds to j/o = and according to (jsj), 
( I55I ) the solution to the CTC becomes singular. From Fig. || we see that apart 
from a small neighborhood around the singular points the CTC gives a good 
description of the 2-soliton train of the MNSE; the singular points match rather 
well with the points at which the two solitons are closest to each other. The 
distance to the first singular points matches T2s^i/4 with T2s,i given by formula 
(H) with yo = 0. 

Choosing the solitons to have different widths leads to 710 7^ in Eq. ( |56| ) 
and removes the singularity of the corresponding solution of the CTC system 
even if F = 0. This can be seen from Fig. || which corresponds to a BSR. Of 
course now the match between the MNSE simulation and the CTC solution is 
better than in the previous case. 

The situation changes if we consider solitons with phase differences such 
that F = TT. There we find a threshold value for di = —d2 = {vi — 1^0) /^o, see 
Eq. (p3|). Whenever di < Acr,2 we get an AFR (see Fig. ||a while for di > Acr^2 



the MNLS 



we get an BSR (see Fig. ||b.) 
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20 40 60 80 100 

Time, t 

Figure 1: Two-soliton interactions and their comparison with the CTC model. 
SoHd curve: numerical results; dashed curves: predictions from the Toda chain 
model, vi = U2 = 1.0, 6i = 0, ^2 = tt. 




Time, t 



Figure 2: Two-soliton interactions and their comparison with the CTC model. 
Solid curve: numerical results; dashed curves: predictions from the Toda chain 
model, vi = 0.95, f2 = 1.05, = 0, ^2 = n. 
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Figure 3: Two-soliton interactions and their comparison with the CTC model. 
Sohd curve: numerical results; dashed curves: predictions from the Toda chain 
model, a) vi ^ 0.97, 1^2 = 1.03, 61 = 0, (S2 = 0; b) = 0.96, 1^2 = 1.04, Si = 0, 
^2 =0. 



Let us now consider the 3-soliton interactions. The choices of the soliton pa- 
rameters illustrates each of the three main configurations outlined in Section 5.2 
above. 

The figures ^ provide examples of 3-soliton configurations with q = char- 
acteristic for case I. Both sets of parameters are such that $ = tt. Besides on 
Fig. ^ we have di < Acr,3 and as a consequence an AFR must follow. In the 
next Fig. ^ we have di > Ad, 3 for which the CTC model predicts a BSR; the 
match with the simulation here is not that good. 

The figures ^ show a 3-soliton configurations with p = characteristic for 
case IL In fig | Id the set of widths is Wl and di < Acr,3 and therefore a MR 
follows 

In the last two figures ^ we used W2 set of soliton widths and a choice of 
parameters characteristic for case III, i.e. p is real while q is purely imaginary. 
In Fig. |a Q > with BSR, and in Fig. |b, we have Q < and AFR. 



7 Conclusions 

In this paper we extend the formalism by Karpman and Solov'ev proposed to 
describe the NSE 2-soliton interaction ||^ and generalized to arbitrary number 
of solitons 1^, ^, 0, js) , to the case of the modified nonlinear Schrodinger equa- 
tion. The aim of our paper was two-fold. First, we would like to investigate a 
possibility to apply an integrable chain-like model to capture adiabatic dynam- 
ics of MNSE solitons within the iV-soliton train. Because a functional form of 
the MNSE soliton is not of the familiar hyperbolic-secant type with a real ar- 
gument, we might expect an existence of some novel features as compared with 
the NSE case. We show that, under specific well-defined conditions, the dynam- 
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Figure 4: Three-soliton interactions and their comparison with the CTC model. 
Solid curve: numerical results; dashed curves: predictions from the Toda chain 
model, a) ui = 1.04, V2 = 1.0, = 0.96, 5i = 0, ^2 = -0.0392, 83 = 0.0016; b) 
z/i = 1.07, U2 = 1.0, V3 = 0.93, ^1 = 0, (52 = -0.0676, 5^ = 0.0049. 




160 

Time, t 



150 

Time, t 



Figure 5: Three-soliton interactions and their comparison with the CTC model. 
Solid curve: numerical results; dashed curves: predictions from the Toda chain 
model, a) 1^1 ^ 1.04, 1^2 = 1.0, 1^3 = 0.96, Si = 0, (52 = 2.1703, S3 = 0.0016; b) 
ui = 1.02, U2 = 0.96, U3 = 1.02, 61 =0,82 = -1.0862, ^3 = 0.0420. 
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Figure 6: Three-soliton interactions and their comparison with the CTC modeL 
Sohd curve: numerical results; dashed curves: predictions from the Toda chain 
model, a) = 1.02, V2 ^ 0.96, V3 = 1.02, Si ^ 0, 82 = 3.142, S3 = 0.0420; b) 
ui = 1.02, 1^2 0.96, 1^3 = 1.02, 61 = 0,62 = 0.0, ^3 = 0.0420. 



ical system of AN equations for soliton parameters is reduced to the completely 
integrable complex Toda chain model with N nodes. This is a strong argument 
in favor of universality of the CTC model for A^-soliton interactions. Though 
the same CTC arises also for the NSE, there are a few peculiarities inherent in 
the MNSE solitons. In particular, we found out more complicated phase be- 
havior of the A^-soliton train. Using the integrability of the CTC, we are able 
to predict various asymptotic regimes of the MNSE 7V-soliton train evolution. 
Besides, we point out the sets of the initial soliton parameters corresponding to 
each of the dynamical regimes. Numerical simulations of the MNSE 2- and 3- 
soliton interactions are in very good agreement with the CTC-based predictions. 
Evidently, the results obtained can be extended to treat also multicomponent 
(vector) generalizations of both the NSE (see, e.g., |3^, Q and references 
therein) and MNSE [||, |||. Work in this direction is now in progress. We 
note that in non-integrable wave systems, Toda-chain type equations may still 
be derived for the adiabatic interaction of N nearly identical solitary waves, but 
such equations are generally non-integrable as well |jl^, ^ . 

Secondly, we consider the MNSE as a true starting integrable model to de- 
scribe subpicosecond pulse evolution in nonlinear media. Strictly speaking, to 
justify a relevance of our results to actual ultrashort pulses, we should also ac- 
count in our model at least two additional effects, the third-order dispersion and 
intrapulse Raman scattering. These effects break the integrability of the CTC, 
and we are faced with a truly perturbed MNSE. Following the lines of recently 
established interrelations between the perturbed NSE and perturbed CTC ||], 
we can extend the above formalism to account for small actual perturbations 
which act along with the effective perturbation (14). The corresponding results 
will be published elsewhere. The single MNSE soliton dynamics in the presence 
of the intrapulse Raman scattering is discussed in the recent paper p9||. 
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Recently we were informed that the CTC model arises also in the case 
of the soliton-train propagation in a system governed by the classical Thirring 
model ^ ||. This seems natural in view of the facts that: i) CTC describes 
the adiabatic soliton interactions for all NLEE of the NLS hierarchy; ii) MTM 
is just another representative of the MNLS hierarchy. 

There remain several natural questions that will be addressed in sequels of 
this paper. The first one is the limit a — > in which we should recover the results 
for the NSE A^-soliton trains. We have proved that the Karpman-Solov'ev-like 
equations for MNSE N solitons transform under this limit to the known NSE 
formulae. The second one concerns the treatment of the perturbed versions 
of the MNSE and the corresponding perturbed CTC model; for the NSE such 
perturbed CTC models have been briefly analyzed in |^ . 
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